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Abstract 

This paper investigates given an assembly graph, find the possible micronuclear sequences 
in terms of MDSs and IESs; which represent Hamiltonian polygonal paths. We will consider the 
O | orientations of the assembly graph and Hamiltonian polygonal path. To obtain a micronuclear 

D sequence of this path, we will compare orientation of the Hamiltonian polygonal path with 

respect to the orientation of the assembly graph. We concentrate on two assembly graphs 
0^ with two 4-valent vertices. These Hamiltonian polygonal paths influence the type of smoothings 

t-H performed. There are two types of smoothings parallel and non-parallel smoothing. We consider 

a micronuclear sequence after one smoothing of a vertex. 

o 

O 1 Introduction 

During the macronuclear development of ciliates (Stylonychia or Oxytricha), DNA recombination 
occurs. Each nucleotide sequence of the macronucleus may be in non-consecutive order in the 
micronucleus as macronuclear destined sequence (MDSs) separated by non-coding DNA. During 
macronuclear development the non-coding fragments (internal eliminated sequence, IESs) that in- 
J> terrupt MDSs in the micronucleus are excised. However, the order of MDSs in the micronucleus 

may not be consecutive, and may be reversed relative to the others pp. The rearrangement of a 
micronuclear sequence can be seen as MDSs and IESs. The assembly graphs represent these mi- 
cronuclear genes. An assembly graph is a finite connected graph where all vertices are rigid vertices 
of valency 1 or 4 |T] . We will focus on assembly graphs that have Hamiltonian polygonal paths that 
represent a single gene combination in the micronucleus, such graphs are called realizable graphs. 
A polygonal path is a path in which every two consecutive edges are neighbors with respect to 
their common vertex or in other words, makes a "90 degree" turn at every rigid vertex pp. A 
. ,— i Hamiltonian polygonal path is a path in which each 4-valent is visited exactly once, makes a "90 

degree turn" at each rigid vertex and no edges are repeated. A transverse path is a path in which 
every two consecutive edges are not neighbors. The given assembly graph (transversal) corresponds 
to the micronuclear DNA segment before recombination. A transverse graph (transversal) is where 
a single path contains every edge is visited and none is repeated, which are called simple assembly 
graphs [3] . We are interested in what types of rearrangements and how many distinct genes can be 
encoded by an assembly graph [4j. 



2 Forming Micronuclear Sequences 

2.1 Notation for Hamiltonian Polygonal Paths and Micronuclear Sequences 

Given a transversal with possible Hamiltonian polygonal paths (HPPs) , we will write possible DNA 
recombinations of micronuclear sequences. In order to understand Hamiltonian polygonal paths 
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and micronuclear sequences, the following will be some notation to help describe them. 

• n represent 4-valent vertices 

• eo, • • • , e2n represent edges, are enumerated by the way they are encountered by the transverse 
path, this fixes an orientation of the edge 

• vq, . . . , v n +i represent vertices of assembly graph 

• r = 71, . . . , 7 n+ i where T is one HPP and 71, 7„+i are ending edges and 72, . . . , j n are full 
edges of G and T, let G assembly graph be n 4-valent vertices 

• P= 72, ■ ■ ■ ,7n 

• Iq, ■ ■ ■ , In+i represent IES sequences, for n is number of vertices in assembly graph 

• Mi, . . . , M n+ i represent MDS sequences, for n is number of vertices in assembly graph 

• Mi, . . . , M n+ i represent inverted MDS sequences, for n is number of vertices in assembly graph 
, dashed lines represent the Hamiltonian polygonal path 

• solid line represent the transversal path 

• A= incoming edges, B= outgoing edges; of transversal path orientation 

We indicate with A an incoming edge incident to a vertex and with B outgoing edge. 

Now, with a given assembly graph and Hamiltonian polygonal path without ending edges (P), 
consider the first and last ending edges of the path, (71,7^+1) and vertices of the path. For P= 
72,. . . ,7 n , which are the complete edges from the assembly graph. Also, HPP= wo 71,. . . , w n -ij n+ i; 
the order of encountering these edges defines an orientation of the HPP. Label the incoming and 
outgoing edges that are neighbors of 71 incident to v « with 1A and IB respectively. Then, label the 
incoming and outgoing edges that are neighbors of 7^+1 incident to V{ n with 2 A and 2B respectively, 
for some i. Each choice of labels 1A, IB or 2A, 2B indicate a possible choice of "start" or "end" 
of an MDS sequence. To denote a specific HPP with chosen ends path, use P(7i,7 n+ i). Each 71 
and 7„ + i will be either one of the labels of ending edges and P is composed of complete edges. 
This path will take upon its own orientation; separate from the orientation of the assembly graph. 
To obtain the micronuclear sequence, compare the orientation of the HPP to the original path 
orientation of the assembly graph. These orientations will indicate the MDSs' orientation. Let, 
A(e;) be the label for a for some i, then the micronuclear sequence corresponds to A(e G ),. . .,A(e2 n )- 
All other edges not visited by HPP are labeled with IESs. Note, an edge may be both 1A and 
2B or IB and 2A (See Figure 7). Such edge is not in the list of HPP edges, 71,. . . ,7n+i- Hence, 
it also has a label Ik for an IES, for some k, so one of the edges may have up to three labels 1A, 
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Ik, 2B for example, see Figure 2. In the micronuclear sequence, if the orientation of 7$ obtained 
from the transversal is opposite, the orientation ^ obtains from HPP, label ji with Mj otherwise, 
its label remains Mi. We assume that the assembly graphs are open ended, so we will always start 
and end with IESs in the micronuclear sequences. But, to find all possible micronuclear sequences, 
we will find all possible combinations of orientations of each HPP with orientations of assembly 
graph. So, each HPP has four distinct combinations of orientations, with respect to orientation 
of the HPP and assembly graph. The four distinct orientations are T, T R , F~ and T~ R . T is an 
HPP, P(7i,7 n+ i). is the opposite orientation of the assembly graph and same HPP orientation. 
r~ is the same orientation as the assembly graph and opposite orientation of HPP. r - ^ is both 
opposite orientations of assembly graph and HPP. To obtain the micronuclear sequences of the 
other 3 distinct orientations, we go through a similar process of the HPPs' of the following figures, 
once the orientations of the assembly graph and HPP are fixed. The following HPPs illustrations 
are T, with its corresponding micronuclear sequence. 

2.2 Possible Micronuclear Sequences Represented by Assembly Graph 1212 




Figure 1: Assembly graph 1212 consists of the two 4-valent vertices and five edges and transverse 
path starts with vo,. . . ,vs. For example, the graph above starts with an edge eo that first visits 
^1)^2)^1)^2, giving us assembly graph 1212. 
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l)e 2 (lA,2A)_ 
I M 1 IiM 2 hM 3 I 3 



2)e 2 {lA,2B) 
I Q MihM 2 hMzh 




Figure 4 Figure 5 



The figures above are HPPs (r) grouped consisting of the path (e 2 ). 

Figure 2 has two vertices so there is 1 full edge in the HPP (e 2 ), with ending edges (1A,2A). 
Full edge, in the sense that it is a complete edge of the assembly graph and T. So, the first 
line of notation above Figure 2 denotes HPP as e 2 (lA, 2 A). The second line of notation above the 
assembly graph represents the micronuclear sequence of the Hamiltonian polygonal path. To obtain 
the micronuclear sequence, read the labels of the transversal path of the given graph. For example, 
we always start with an IES (Iq) following 1A, which is the first MDS (Mi) and continues along 
ei, which is I\. Now, moving along e 2 (M 2 ) we continue to I 2 which is part of e 3 . Then, with M3 
(2 A), however, we must complete the assembly graph and micronuclear sequence with an IES (Is). 
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Notice, the MDS sequences in the HPP are inverted, because their orientations are the opposite 
with respect to the transversal. Obtaining the micronuclear sequences for T R , T~ and r - ^ is the 
same process as obtaining T, as long as the combinations of orientations for each are distinct, with 
respect to the assembly graph and HPP. 

5)e 3 (lA,2A) 6)e 3 (lA,2fl) 
I M 1 I 1 M 3 I 2 M 2 I 3 I M 1 I 1 M 3 I 2 M 2 h 




M 2 



Figure 8 Figure 9 

The figures above are HPPs (r) grouped consisting of the path (es). 
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9)ei(L4,2,4) 
/0M2/1M1/2M3/3 



10)ei(lA,25) 
J0M2I1M1J2M3I3 



M. 



1A 



M. 



-M 1 



1B 



'2A 



M, 



1A 



1B 



2B 



Figure 10 



Figure 11 



ll)ei(LB,2A) 
I M2liMihMh 



12)ei(LB,2B) 
I M 2 IiM 1 I 2 M 3 





Figure 12 



Figure 13 



The figures above are HPPs (r)grouped consisting of the path (ei} 
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HPP 


Micronuclear Sequence 


Inverted Micronuclear Sequence 


1) e 2 (lA,2A) 

2) e 2 (lA,2B) 

3) e 2 {lB,2B) 

4) e 2 {lB,2A) 


I G MihM 2 I 2 M 3 I 3 
/0M1/17W2/2M3/3 
/0M1/1M2/2M3/3 
hWiIi~M 2 ~I 2 W 3 I 3 


I0M3J1M2/2M1J3 
/0M3/1M2/2M1/3 
/0M3/1M2/2M1/3 
I0M3J1M2/2M1J3 


B)e 3 {lA,2A) 

6) e 3 (lA,2B) 

7) e 3 (lB,2A) 

8) e 3 (lB,2B) 


I MihM 3 I 2 M 2 I 3 

/0M1/1M3/2M2/3 
I MihM 3 I 2 M 2 I 3 
IqW 1 I 1 M 3 I 2 M 2 I 3 


IoM^MshMih 
/0M2/1M3/2M1/3 
I M 2 IiM 3 I 2 M 1 I 3 


9) ei(lA,2A) 

10) ei(lA,2S) 

11) ei(lB 2A) 

12) ei(15 2B) 


/0M2/1M1/2M3/3 
/0M2/1M1/2M3/3 

hM 2 hMiI 2 M 3 I 3 
/0M2/1M1/2M3/3 


/0M3/1M1/2M2/3 
/0M3/1M1/2M2/3 
hM 3 hMiI 2 M 2 I 3 



The first column in the table above are the possible HPPs ( T) of assembly graph 1212. The 
second column represent the micronuclear sequences obtained from the HPPs. The third column 
is r^. Some sequences in T~ and r - ^ corresponded to T, T^, which leaves us with twenty four 
distinct labellings of the assembly graph 1212. Notice, that T R is complementary to the original, 
r and therefore represent the same DNA sequence [2j. Graphically, the difference between T R and 
r is the orientation of the assembly graph, which inverts the sequences when obtained. 
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2.3 Possible Micronuclear Sequences Represented by Assembly Graph 1221 




Figure 2: The assembly graph above is 1221, which is formed similarly to the previous assembly 
graph. 



l)ei(lA2i) 
I M 2 hMzI 2 M l h 




Figure 2 



2)ei(lA,2B) 
I M 2 hM 3 I 2 M 1 I 3 




Figure 3 



3)ei(lB,2i4) A)e l {lB,2B) 
I0M2I1M3I2MI3 70M2J1M3J2M1/3 




Figure 4 Figure 5 



The figures above are HPPs (r) grouped consisting of the path (ei). The corresponding mi- 
cronuclear sequences are obtained in a similar fashion as the previous micronuclear sequences of 
previous assembly graph 1212. 
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5)e 2 (lA,2B) 

I Q MihMzI 2 M2h 




Figure 6 



6)e 2 (lA_2,4) 

/0M1/1M3/2M2/3 




Figure 7 



hM x hM z I 2 M 2 h 




8)e2(_lB, 2B) 

I Q M 1 I 1 M 3 I 2 M2h 




Figure 8 



Figure 9 



The figures above are HPPs (r) grouped consisting of the path (e 2 ). The corresponding mi- 
cronuclear sequences are obtained in a similar fashion as the previous micronuclear sequences of 
previous assembly graph 1212. 
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HPP 


Micronuclear Sequence 


Inverted Micronuclear Sequence 


1) ei(L4,2A) 

2) ei(lA,2B) 

3) ei(15,2A) 

4) 61(15,25) 


/0M2/1M3/2M1/3 
I§M 2 hM 3 I 2 M x I 3 
/0M2/1M3/2M1 
IqM 2 I x M 3 I 2 M x 


/0M1/1M3/2M2/3 
/0M1/1M3/2M2/3 
M1/0M3/1M2/2 
M1/0M3/1M2/2 


5) e 2 (lA,2B) 

6) e 2 (lA,2A) 

7) e 2 (lB,2A) 

8) 62(15,25) 


I M 1 I 1 M 3 I 2 M 2 I 3 
IoMilJvhhM^h 

IoWihMzhlvhh 


I M 2 hM 3 I 2 M 1 I 3 

I M 2 hM 3 I 2 Mh 
hM 2 hM 3 I 2 M x l 3 
I M 2 I 1 M 3 ~I 2 M 1 I 3 



The first column in the table above are the possible HPPs ( T) of assembly graph 1221. The 
second column represent the micronuclear sequences obtained from the HPPs. The third column 
is T R . Notice, the sequences of T R are complementary to the DNA sequences of T. Which leaves 
eight distinct micronuclear sequences of this fixed orientation. For example: 
The sequence of line 8 is the sequence in line 1 
The F R sequence of line 7 is the sequence in line 2 
The F R sequence of line 5 is the sequence in line 3 
The F R sequence of line 6 is the sequence in line 4 

But, we have to find all possible micronuclear sequences for assembly graph 1221, which also 
considers T~ R and T~. The other eight distinct micronuclear sequences come from T~ R and T~ as 
follows, 



HPP Orientation 


Micronuclear Sequence 


r- 


hM 2 l x M x I 2 M 3 I 3 
I M 2 hM 1 I 2 M 3 I 3 
I M 2 ~hM 1 I 2 M 3 I 3 
1^M~ 2 1{M{1 2 M 3 1 3 


r-ji 


hM 3 I x M x l 2 M 2 I 3 
IqM 3 I x M x 1 2 M 2 I 3 
IqM 3 I x M x 1 2 M 2 1 3 
IqM 3 I x M x 1 2 M 2 1 3 



with a total of sixteen distinct micronuclear sequences of assembly graph 1221. 
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3 Smoothing of Assembly Graphs 



Smoothings of assembly graphs represent alignment of pointers (MDS pointers) in the micronuclear 
sequences and the rearrangement. Each vertex, there are edges incident to the vertex labeled with 
MDSs and IESs that cross each vertex. When smoothed at the vertex (removing the vertex), the 
pointers align (See Figure 10). During this process, MDSs align and the in the micronuclear se- 
quence, recombine MDSs and IESs. There are two types of smoothing, N (non-parallel) smoothing 
and P (parallel) smoothing (See Figure 11). If the pointer sequences are aligned in parallel, then 
the smoothing at the vertex follows the predetermined direction of the assembly graph (parallel 
smoothing). If the pointer alignment is anti parallel, then the smoothing of the vertex changes 
the path of orientation of the original graph [lj. Smoothing of each vertex in an assembly graph, 
is determined by the HPPs of the assembly graph. The type of smoothing is predetermined by 
the given HPP, in such a way the HPP stays connected. The writing of the "new" micronuclear 
sequences after the smoothing, is similar to writing of the micronuclear sequence before the smooth- 
ing. Use similar notation from Sec. 3 to write "new" micronuclear sequences. The difference from 
the original transversal " old" micronuclear sequence to the new assembly graph " new" micronuclear 
sequence (not necessary a transversal) is recombined MDSs, since that's what smoothing does to 
the sequences. Notice, the appearance during P smoothings, the assembly graph is broken up into 
components, that is why there are bracket notations. The brackets show a separation in MDSs 
within the assembly graph. 




p-smoothing n-smoothing 

Figure 11 
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3.1 Smoothing of Assembly Graph 1212 



The following figures represent smoothings of HPPs on assembly graph 1212. Where corresponding 
micronuclear sequences of the HPPs is on left side of arrow. The smoothings now correspond 
to a new micronuclear sequence, which is represented on the right side of the arrow. As stated 
previously, the HPPs predetermines the type of smoothing of the HPPs on the assembly graph, 
which is one reason why the sequences are grouped in this manner. The other is with respect to 
which vertex is removed. Notice, the merging of MDSs after vertex removal. 



N smoothing at vertex 1 



1) J0M1J1M2J2M3I3 

2) /0M1/1M2/2M3/3 
7) /0M1/1M3J2M2/3 
8) 

11) I M 2 I 1 M 1 I 2 M 3 I 3 

12) I M 2 I 1 M 1 I 2 M 3 



; Mi )2 / 1 M 3 / 2 
I M 1>2 hM 3 I 2 

IoM3,i,2h 
I M 3 hM h2 I 2 

loM 2 ,i, 3 h 

I M 2>1 hM 3 I 2 




Figure 3: The following smoothing corre- 
sponds to the first micronuclear sequence of 
assembly graph 1212. 



P smoothing at vertex 1 



3) IoMihAhhMs 

4) I Mi 1 1 M^I 2 M 3 1 3 

5) I MihM 3 I 2 M 2 I 3 

6) hMxhM^hNhh 

9) hM 2 I x M x I 2 M 3 I 3 

10) hM 2 I x M x I 2 M 3 I 3 



/0M3, [hM 2 ,ih] 

I M 3 h, [M^ih] 

I M lj2 h[M 3 I 2 ] 

IoMi^hlM^] 

I M 3 h[M lj2 I 2 ] 

I Q M 3 [hM lt2 I 2 } 




Figure 4: The following smoothing corre- 
sponds to the third micronuclear sequence 
of assembly graph 1212. 



As seen in Figure 3, represents the assembly graph for first micronuclear sequence of assembly 
1212, when vertex one is removed. We can see in the micronuclear sequence before the vertex 
removal, had separate MDSs. After, the vertex removal, there were merging of MDSs, Mi and 
M 2 to M\ y2 - To obtain the micronuclear sequence, read the transversal path with respect to the 
smoothed HPP, considering the fixed orientations of both. 
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P smoothing at vertex 2 



N smoothing at vertex 2 



1) /0M1/1M2/2M3/3 

4) /0M1/1M2/2M3/3 

6) I MihM 3 I 2 M 2 I 3 

8) /0M1/1M3/2M2/3 

10) I Q M 2 hMxI 2 Mzh 

12) I Q M 2 I 1 M 1 I 2 M 3 I 3 



7 Mi/i[%/ 2 ] 

J MlIl[M3,2/2] 

I Mi/i[M 3)2 J 2 ] 
I M 2)3 [/iMiJ 2 j 
/ M 2)3 [7iM 1 / 2 ] 



2) I MihM2hM 3 I 3 

3) hM x hM2l 2 M 3 I 3 
5) I MihMzhM 2 I 3 

7) /0M1/1M3/2W3 

9) I M 2 I 1 M 1 I 2 M 3 I 3 

11) hM 2 l x M x l 2 M 3 l 3 



I MihM 2 , 3 I 2 
IpMxhM^h 
I MihM 3>2 I 2 
IqM 1A2 I]_ 

/ M2,3/lMiJ 2 

/oAf 2l3 ,i/i 



• M 2,3 



M, 




Figure 5: The following smoothing corre- 
sponds to the first micronuclear sequence of 
assembly graph 1212. 



Figure 6: The following smoothing corre- 
sponds to the second micronuclear sequence 
of assembly graph 1212. 



As seen in Figure 5, represents the assembly graph for first micronuclear sequence of assembly 
1212, when vertex two is removed. We can see in the micronuclear sequence before the vertex 
removal, had separate MDSs. After, the vertex removal, there were merging of MDSs, M 2 and M 3 
to M2,3. Notice, there are 2 components of the transversal. To obtain the micronuclear sequence, 
read the transversal path with respect to the smoothed HPP, considering the fixed orientations of 
both. The brackets are used to indicate the beginning of the second component of the transversal, 
with respect to the "new" micronuclear sequence. 



3.2 Smoothing of Assembly Graph 1221 

The following figures represent smoothings of HPPs on assembly graph 1221. Where correspond- 
ing micronuclear sequences of the HPPs is on left side of arrow. The smoothings now correspond 
to a new micronuclear sequence, which is represented on the right side of the arrow. As stated 
previously, the HPPs predetermines the type of smoothing of the HPPs on the assembly graph, 
which is one reason why the sequences are grouped in this manner. The other is with respect to 
which vertex is removed. Notice, the merging of MDSs after vertex removal. Obtaining the "new" 
micronuclear sequences is similar to smoothings of assembly graph 1212. 
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P smoothing at vertex 1 



N smoothing at vertex 1 



1) I M 2 IiM 3 I 2 M 1 I 3 

2) IoAhhMshMih 

7) /0M1/1M3/2M2/3 

8) Iq Mi 1 1 M 3 1 2 M 2 I 3 



I M h2 hM 3 I 2 
I M 1>2 h M 3 I 2 

hM 2X3 h 
IoM 2>1 hM 3 I 2 



3) I M 2 IiM 3 I 2 MiI 3 

4) IoM 2 hM 3 I 2 MiI Q 

5) /0M1I1M3/2M2/3 

6) IoM-LliMshM-zh 



I M 3 hM 2A I 2 
IoMshM^xh 
I M 1>2 I 1 M 3 I 2 
I M 1<2 I 1 M 3 I 2 



M , 




Figure 7: The following smoothing corresponds 
to the first micronuclear sequence of assembly 
graph 1221. 




Figure 8: The following smoothing corre- 
sponds to the third micronuclear sequence 
of assembly graph 1221. Removing vertex 1 
allowed M\ and M 2 to align forming M 2 , 1. 



P smoothing at vertex 2 



2) I M 2 hM 3 I 2 MiI 3 

4) I M 2 hl^I 2 MiI 3 

6) I0M1/1M3/2M2/3 

7) lQM x hM 3 l 2 M 2 l 3 



I M 2>3 hMxI 2 
IpM^hMxh 
J M 1 J 1 M 3i2 / 2 
hMxhM^h 




Figure 9: The following smoothing corresponds 
to the second micronuclear sequence of assembly 
graph 1221. 



N smoothing at vertex 2 



1) hM 2 hM 3 l 2 MxI 3 

3) I M 2 hM 3 I 2 MiI 3 

5) hMihM^M^h 

8) IoMJiMshM-zIs 



I Q M 2t3 hMiI 2 
I Q M 2 , 3 IxMxI 2 
IpMihM^h 
I M 1 I 1 M 3t2 I 2 




Figure 10: The following smoothing corresponds 
to the second micronuclear sequence of assembly 
graph 1221. 



4 Conclusion 

The forming of micronuclear sequences through HPP can influence the resulting sequences. Such 
as, the smoothing of vertices, which behave in a manner to preserve the HPP. But, at the same 
time eliminate some IESs and align MDSs. The previous figures were for two vertices, as the 
vertices increases, the HPP gets more complicated. The micronuclear sequences are longer, due 
to the fact more edges are used to form the HPP. The smoothings are more abstract; since the 
assembly graphs are larger and more smoothings since there are more vertices. There must be 
four distinct micronuclear sequences for a given path P(P(1A, 2A), ...). Also, given an orientation 
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of the transversal, each sequence has an inverse. Notice, the second table on page 10 did not 
include inverses, since they corresponded with the micronuclear sequences of the HPP orientations. 
As an observation, there are at least eight micronuclear sequence (unless they are complementary 
or assembly graph 11). For the other orientations F~ and r - ^ there are eight new sequences, 
potentially. 
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